function g = gradf1(x)
    global G b;
    g = G * x + b;   
end